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ABSTRACT 

We present pulse-phase X-ray spectroscopy of the high-mass X-ray pulsar SMC X-1 from five differ- 
ent epochs, using Chandra X-ray Observatory ACIS-S and XMM-Newton EPIC-pn data. The X-ray 
spectrum consistently shows two distinct components, a hard power law and a soft blackbody with 
kT^B ~ 0.18 keV. For all five epochs the hard component shows a simple double-peaked pulse shape, 
and also a variation in the power law slope, which becomes harder at maximum flux and softer at 
minimum flux. For the soft component, the pulse profile changes between epochs, in both shape and 
phase relative to the power law pulses. The soft component is likely produced by reprocessing of the 
hard X-ray pulsar beam by the inner accretion disk. We use a model of a twisted inner disk, illumi- 
nated by the rotating X-ray pulsar beam, to simulate pulsations in the soft component due to this 
reprocessing. We find that for some disk and beam geometries, precession of an illuminated accretion 
disk can roughly reproduce the observed long-term changes in the soft pulse profiles. 
Subject headings: accretion, accretion disks — stars: neutron — X-rays: binaries: individual (SMC 
X-1) — stars: pulsars: individual (SMC X-1) 



1. INTRODUCTION 

The general picture of accretion onto X-ray pulsars, 
consisting of a flow in a disk to the magnetosphere and 
along the dipole fleld to the magnetic poles of the neutron 
star, has been known for decades (sec Nagasc 1989, and 
references therein). However, despite extensive study, 
some important aspects of this process remain myste- 
rious. In particular, the details of accretion near the 
magnetosphere, where the neutron star's magnetic field 
begins to do minate the flow, i s a complex and challeng- 
ing problem llElsner fc LambI 119771 iGhosh et all Il977t 
IChosh fc T.a,mblIT97i: 

1.1. Superorbital variation in X-ray pulsars 

One way to explore an accretion flow observationally 
is with systems where the shadow of the disk causes the 
X-ray source to vary with time. There are now ~20 
known accreting binary systems that exhibit superor- 
bital periods, some of which are caused by shadowing of 
the s ource flux by an accretion disk (see Clarkson et alj 
120031 and references therein). Of these systems, a small 
group of X-ray pulsars (SMC X-1, LMC X-4 and Her- 
cules X-1) show luminous, pulsing X-ray emission which 
allows for study of the accretion flow in unusual de- 
tail. These systems each show a 30-60 d superor- 
bital variatio n, attributed to a warp e d, precessing ac- 
cretion disk llTananbaum et alJ Il972t iJones &: FormanI 
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In the low- mass system Her X-1, the precessing disk 
has been studied extensively. Regular 35-day varia- 
tions have been observed in the X-ray flux, pulse pro- 
files, and spectrurn, as well as optical and UV features 
fe.g.. iVrtilek et all 11994 iDeeter et af] 119981 and refer- 
ences therein). In particular, occulation by the disk pro- 
duces variations in the pulse profile shapes, which can be 
used to infer the overa ll warped disk geometry ijLeahvl 
l200llScott eralll200(il) . 



SMC X-1 and LMC X-4 differ from Her X-1 in that 
they have high-mass (O and B star) companions, hence 
we expect the accretion mechanism to be somewhat dif- 
ferent due to the presence of powerful stellar winds. How- 
ever, the high X-ray luminosities of SMC X-1 and LMC 
X-4 (> 10^^ ergs s""'^) require Roche lobe overflow in ad- 
dition to accretion from the wind, and there is evidence 
for persistent, warped accretion disks. 

For example, in LMC X-4, precession of the disk is 
believed to cause changes in the X-ray power-law slope 
and flouresccnt iron line flux as a function of superorbital 
phase (Naik & Paul 2003, 2002) , as well as v ariations in 
the ultraviolet spectra (Prcc iado et a l. 2002). For SMC 
X-1, observations of X-ray spectra ( Woidowski et all 
119981) and ohotoionized lines in the UV and X-ray 
() Vrtilek et al.l 120011 12005|) give evidence for shadowing 
of the central neutron star by the precessing d isk. From 
the observed line features. IVrtilek et a l. (2005) conclude 
that the occulation of the source is caused by the inner 
regions of the disk, near the magnetospheric radius ^10^ 



1.2. Disk reprocessing and the soft X-ray emission 

In Her X-1, reprocessing of the hard X-rays from the 
neutron star by the inner accretion disk produces a soft 
(< 1 kcV) X-ray component that pulses out of phase 
with the ha rd (> 2 keV) pulses (e.e. lEndo et al.li2000t 
iRamsav et al. 2002). Thus the X-ray beam from the 
neutron star acts as a "flashlight" that illuminates the 
inner disk (Fig. 1). 

Many luminous X-ray pulsars, including LMC X-4 and 
SMC X-1, have similar soft spectral com ponents that 
are li kely also caused by disk reprocessing (jHickox et alJ 
120041) . We expect the flux of such a soft component to 
pulsate as the pulsar beam sweeps around and illumi- 
nates the disk. The soft pulses will have the same period 
as those from the hard X-ray beam, but will have a dif- 
ferent shape and relative phase. Moreover, precession 
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Fig. 1. — Schematic of disk reprocessing in X-ray pulsars. The 
hard (power law spectrum) X-ray beam sweeps around and illumi- 
nates the accretion disk, which re-radiates a soft X-ray (blackbody) 
component. The observer sees pulses from both the power law and 
blackbody components, but these differ in pulse shape and phase. 



of the accretion disk will change the observer's view of 
the inner region, and so the shape and phase of the re- 
processed component should vary with superorbital pe- 
rio d. Hints of such an effect have been seen in Her X- 
1: iZane et all l|2n04ll showed using XMM data that the 
phase offset between the soft (0.3-0.7 keV) and hard (3- 
10 keV) pulse profiles vary with time, and appears to 
be correlated with precession of the disk. A soft compo- 
nent that pulsates separately f rom the hard compo nent 
has also be en seen in LMC X-4 llNaik fc Paul2004bD and 
SMC X-1 llPaul et al.ll20(11 INaik fc Pau]ll2004aH . 

1.3. SMC X-1 

SMC X-1 was first identified as a p oint source in Uhuru 
observation s bvlLeong et all il97lD . Eclipses were dis- 
co vered bv ISchreieretld.l l)1972() ~and X-ray pulsations 
by iLucke et^dT " i[T9'7fiD . The optical counterpart is S k 
160, a BO I supergiant l|Webster et al...l97 2: Liller .,1973t) . 
The orbit is close to circular, with orbital period ^^3.9 
days and pulse period 0.7 s. The source has a high X-ray 
luminosity, Lx ~ 3 x 10^® ergs s~^, in the bright state. 

In SMC X-1 the perio d of the superorbital cyc le oscil- 
lates from 40 to 60 days JWoid owski et al.lll998j) . a vari- 
ation that may itself be periodic on the order of years 
ijClarkson et al.ll2'003D . The high-state, out-of-eclipse X- 
ray spectrum is usually fitted as a cutoff power law with 
r ^ 0.9, an iron line at ^6.4 keV, and a soft component 
modeled as a broken power law, thermal bremsstrahlung, 
or a blackbody with fcTeB 0.16 keV. lPaul et"an 1)20021) 
fit four different models to the soft component and found 
that all gave satisfactory fits. They found that the 
soft component flux pulses roughly sinusoidally, out of 
phase with th e sharper-p eaked pulses from the power 
law. iNcilscn e t al.l l|2004f l created hard and soft pulse 
profiles for ROSAT, Chandra, and XMM data of SMC X- 
1, and found that for different epochs the the soft pulse 
shapes vary while the hard pulses remain similar. 

In this paper we seek to explore the variation in these 
spectral components further, by performing pulse-phase 
spectroscopy for SMC X-1 at several epochs. We have 
also constructed a model to simulate hard and soft pulse 
profiles that would be produced by an X-ray pulsar beam 
illuminating a warped, precessing accretion disk, and 
compare the results with observations. 
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Fig. 2.— RXTE ASM lightcurve for SMC X-1 with arrows 
designating several observations. Quick-look results as provided 
by the RXTE/ ASM team. 
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Fig. 3.— XMM EPIC-pn light curve (0.15-12 keV) for observa- 
tion XlOl, showing entrance to eclipse. Times with fiux > 41 cts 
s~^ (shaded) were included in the spectral analysis. 



2. OBSERVATIONS AND DATA REDUCTION 

We have used three Chandra and two XMM observa- 
tions of SMC X-1, taken in the high state of the superor- 
bital cycle. Details of the observations are given in Table 
1. Pulse prof ile analysis for these data are presented in 
INeilsen et al.l (|2004) , and we will use their abbreviations 
to denote the observations (C102, XlOl, etc.). Fig. 2 
shows the observations on the long-term light curve of 
SMC X-1 from the All Sky Monitor (ASM) on RXTE. 
From this we estimate the superorbital phases given in 
Table 1, taking (pso = at the start of the superor- 
bital high state. These estimates have uncertainty of ± 
a few days or 0go ^ 0.05, because the superorbital in- 
tensity variation in SMC X-1 is quasi-periodic in nature 
and hence does not have a well-defined ephemeris. 

The Chandra observations were each '--^6 ks in duration, 
and they were taken with the ACIS-S3 chip in Continu- 
ous Clocking (CC) mode, in which charge is continually 
read out down the chip. This allows for 3.5 ms time reso- 
lution, but one dimension of spatial information is lost so 
that the X-ray image is simply a straight line. There are 
no other bright sources within the ACIS field of view near 
SMC X-1, so we can determine source and background 
regions from the 1-D image and extract events normally. 
For each observation the source region was a rectangle 
of size 0'.'49 x 1'.'97 centered on the source. Background 
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TABLE 1 

Chandra AND XMM OBSERVATIONS OF SMC X-1 



Observation ID 
(Reference) 


Mission 


Start Time 
(JD-2400000) 


Total 
Exp. (ks)^ 


<f>o-i:h° 


4>so 


Fluxes (10^ 
Power law 


" ergs cm s 
Blackbody 


.1)C 

Total 


400102 {C102) 


Chandra 


51832.17 


6.2 


0.19 


0.16 


9.3 


0.9 


10.2 


400103 {C103) 


Chandra 


51833.34 


6.1 


0.49 


0.16 


9.1 


1.6 


10.7 


400104 {C104) 


Chandra 


51834.31 


6.5 


0.74 


0.16 


9.2 


1.4 


10.6 


0011450101 (XlOl) 


XMM 


52060.59 


27.9 


0.89 


0.37 


6.5 


0.4 


6.9 


0011450201 (X201) 


XMM 


52229.65 


1.6 


0.34 


0.42 


10.1 


1.3 


11.4 



^Clcan, out-of-cclipse exposure. 
^Orbital phase at start of exposure. 
*^Unabsorbed fluxes in the range 0.5—10 keV. 



events came from two rectangles 0'.'49 x 9'.'84 on either 
side of the source. 

For the XMM observations we have used data from 
the EPIC-pn CCD camera, which has the highest time 
resolution of the EPIC cameras. X201 had a clean ex- 
posure of ~27 ks in duration, and was taken in Small 
Window mode, which has 6 ms time resolution. For this 
observation we have extracted source events from a circle 
of radius 51'.'2. XlOl was taken in Full Window mode, 
which has 73.5 ms time resolution. While the total ex- 
posure for XlOl was 40 ks, the bulk of this observation 
was during full e clipse, with only the f irst few ks showing 
substantial flux. iNeilsen et al.l ()2004D analyzed the first 
4.9 ks of data, which consisted of the entire entrance 
into eclipse. However, they found that a partial covering 
absorption model was needed to explain the X-ray spec- 
trum, perhaps due to occultation by the outer layers of 
the star. To avoid this complication, we restricted the 
present analysis to the 1.6 ks of data for which the 0.15- 
12 keV count rate was > 41 counts s~^ (see Fig. 3). The 
background surface brightness for the XMM data was 
< 0.1% of the source brightness, so we did not perform 
background subtraction. 

3. DATA ANALYSIS 

3.1. Dealing with event pileup 

In high count-rate data such as these, the X-rays mea- 
sured by CCD detectors are affected by pileup, caused 
when two or more photons are coincident on one event- 
detection cell during one time-resolution element, or 
frame time fOavis 2001). In this case the detector can- 
not distinguish the multiple events and instead detects 
a single event with a pulse height roughly equal to the 
sum of the pulse heights from the individual photons. 
Pileup has several effects on the detected X-ray events, 
most importantly that the overall detected count rate is 
lowered, and the spectrum is shifted from lower to higher 
energies, due to the absence of real lower energy events 
and the addition of "fake" higher-energy events. 

There is also an effect on the "grade" of the detected 
events. The grade of an X-ray event is a number that 
refers to how the charge cloud from a photon hitting a 
CCD pixel is spread across neighboring pixels. Real X- 
ray events have characteristic grades that can be used to 
distinguish them from unwanted events such as particle 
background interactions. However, when two or more 
photons are piled up, they often result in grades that 
do not match those of individual photons, and so are 
thrown away by this filtering, thus decreasing the count 



rate. Taking into account this "grade migration" is an 
important part of modeli ng pileup. F or a discussion of 
these effects in detail, see lDavi3 (l200l . 

In this study we deal with pileup in three different 
ways: short frame times, excluding events from the cen- 
ter of the PSF (for XMM), and pileup modchng (for 
Chandra). All five observations have short frame times 
<75 ms, which mitigate pileup by reducing the likelihood 
that multiple photons will be coincident in one time res- 
olution element. 

For the XMM data, we tested for pileup using the SAS 
tool epatplot, which compares the distribution of events 
of different grades with that expected for no pileup. 
When events from the full PSF were included, these dis- 
tributions differed significantly, indicating pileup. How- 
ever, because count rates (and thus pileup effects) are 
greatest at the center of the PSF, we found that by ex- 
cluding events from the central 1 /3 radius of the source 
region, we could eliminate any significant pileup effects. 

For the Chandra data we could not simply exclude the 
inner extraction region. Because the PSF is smaller than 
for the XMM and because the use of CC mode concen- 
trates all the events along one dimension, most source 
events are concentrated in a few pixels. Therefore, we 
fit the data in the ISIS package using the pileup kernel, 
which models and corrects for pileup effects in spectral 
fitting (Davis 2001). This kernel includes a grade migra- 
tion parameter a which can be allowed to float to best 
characterize the effect that pileup has on the grade selec- 
tion of the observation. On average 17% of the photons 
in the ACIS spectra arc piled up. This level of pileup 
can be corrected using the ISIS kernel while measuring 
the spectral parameters. 

3.2. Pulse-phase spectroscopy 

Pulse phase-resolved spectra were extracted in the 
range 0.6-10 keV. For the Chandra data we excluded 
the energy range 1.7-2.8 keV in our spectral fits (as in 
INeilsen et al.l 120 04') . to avoid instrumental effects near 
the Au and Si edges. For both XMM and Chandra, we 
cleaned the data for flaring episodes before spectral ex- 
traction. 

The best fit spectral models for the Chandra and XMM 
observations consisted of a power law (P ^ 0.9) plus 
blackbody A:Tbb ~ 0.18, both modified by neutral ab- 
sorption. We also included a high-energy exponential 
cutoff on the power law, as has been seen in Cinga and 
ASCA observations ijWoo et al.lll995t IpTuI et al...2002.'l . 
X201 was the only observation with sufficient high-energy 
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Fig. 4. — Pulse phase-resolved spectra for SMC X-1, for (a) C103 
(b) XlOl, (c) X201. For each observation two spectra are shown, 
from pulse phases that show substantial differences in the spectral 
parameters. The individual blackbody (BB) and power law (PL) 
components are also shown. In the Chandra data, these individual 
components do not add to the total, because of the effect of the 
pileup model. Also note that the energy range 1.7-2.8 keV was 
excluded from the Chandra fits. 



counts to constrain this cutoff, so for the other obser- 
vations we fixed the parameters to the phase-averaged 
value for X201 (£;cut = 6.1 keV and £:foid = 6.8 keV). 
The models can be written as: 



Chandra: f{E) 



[fBB{E) + fpL{E)U^t{E)] 



XMM: J{E) = e 



^ p-cy{E)Nii 



[fBB{E)+fpL{E)f,ut{E)] 



where /bb, /pl, /cut, are the standard XSPEC black- 
body, power law, and high-E' cutoff models. A'piic is the 
ISIS pileup kernel. 

In general, for phase-a veraged spectral par ameters we 
use the results of fits by INeilsen et al.l (|2004') . However 
for XlOl we performed a new phase-averaged fit, because 
we have revised the time interval over which the spectrum 



was extracted (see § 2). We fixed the high-energy cutoff 
as above, giving Nh = 1.3 x lO^^ cm'^, /cTbb = 0.17 
keV, and T = 0.79, with XSPEC normalizations Kbb = 
8.8 X lO^"* (in units of 8.3 x 10"^ ergs cm'^ s'^) and 
JCpL = 3.5 X 10-2 (photons keY'^ cm^^ s-1 at 1 keV). 

To extract phase-resolved spectra, we first corrected 
the photon arrival times for the motion of the Earth and 
the observatory, as well as the orbital mo tion of S MC X- 
1 using the ephemeris found by FWo i dowsld_g^_^ (Q_998). 
Using the pulse periods found bv INeilsen et al.l l|2004l|') . 
we extracted spectra in 10 bins of (/)puise for each of the 
observations. We fit the spectra in ISIS for Chandra and 
XSPEC for XMM, using the models given above. For the 
Chandra fits, we fixed the grade migration parameter a 
at its phase-averaged value. 

For each fit we set the initial parameters to the phase- 
averaged values, and fixed the A^h and high energy cut- 
off while allowing the other parameters to vary. For 
XlOl, the limited statistics required that we also fix the 
blackbody temperature, fcTBB, in order to measure the 
variation in the soft component. Sample phase-resolved 
spectra for C103, XlOl, and X201 are shown in Fig. 4, 
for which changes in shape between spectra at different 
0puisc are clearly visible. The phase-resolved fits have 
xl values of 0.8-1.0 for Chandra, 1.0-1.3 for XlOl, and 
0.9-1.2 for X201. 

We found that the best-fit spectral parameters vary 
significantly with (/ipuise- For each phase we calculated 
the 0.5-10 keV unabsorbed flux in both the blackbody 
and (cutoff) power law components. The 0puiso variation 
of the fluxes, along with F and /cTbb, are shown in Fig. 
5. We flnd that: 

1. At every epoch, the power law flux shows the famil- 
iar double peaked profile, similar to the hard pulse 
profile observed many times previously. The phase 
offset and relative intensities of the two peaks vary 
slightly between epochs. 

2. The power law photon index varies significantly, 
inversely with the power law flux. The spectrum 
gets harder as the power law flux increases, and 
softer as the flux decreases. 

3. For all observations, the blackbody flux shows a 
single main peak, and in the Chandra observations 
there is a shoulder after the main peak. 

4. The shapes and relative phases of the blackbody 
pulses change with time. However the three Chan- 
dra observations (at (/lorb = 0.19, 0.49, and 0.74) 
have very similar soft pulses. This implies that 
the changes in the soft profiles are not correlated 
with orbital motion of the source, and occur on 
timescales longer than a few days. 

To test whether the observed variation is dependent 
on the particular choice of soft component model, we 
repeated the pulse-phase spectroscopy using a thermal 
bremsstrahlung (TB) model in place of the blackbody. 
This fits are acceptable for kTrs ^ 0.3 keV, with very 
similar to the blackbody fits. This is not a phys- 
ical model, because a large, diffuse thermal gas cloud 
could not produce the lumin ous, pulsing soft emission 
observed ijHickox et al.ll200^ . but it allows us to check 
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Fig. 5. — Pulse phase variation of spectral parameters for SMC X- 
1, for (a) Chandra: solid: C102, dashed: C103, dot-dashed: C104 
(b) XlOl, (c) X201. Shown are the relative fluxes from the power 
law and blackbody components, the power law index F, and the 
blackbody temperature A;Tbb in keV (note that the temperature 
was fixed in the XlOl fits). Two pulses are shown for clarity. Errors 
are la, and typical error bars for Chandra are given to the left of 
the plots. The average fluxes in the power law and blackbody 
components are given in Table 1. 



that the measurement of variation in the soft component 
is robust. The soft and hard component flux profiles for 
this model for one observation (C103) are given in Fig. 
6, and are compared to those for the blackbody model. 
The shapes of the profiles differ by < 0.5a, for the two 
models, except for the soft profile in C102, for which they 
differ by ~ 1.2a. Since these differences are not highly 
significant, we conclude that the hard and soft flux pro- 
files are essentially independent of the model chosen for 
the soft component. 

4. PHOTON INDEX VARIATION 




Pulse phase 

Fig. 6. — Pulse-phase variation in spectral parameters for C103, 
using a blackbody (histogram) and a thermal bremsstrahlung 
(dashed line) soft component model. Note that the hard and soft 
flux profiles do not differ significantly for the two models. 



The observed changes in the power law index with 
pulse phase have not bee n previously repo rted for SMC 
X-1; the ASCA study of iPaul et all l(2002.1 found that 
r ~ 1.0 for four different pulse phases. We have in- 
vestigated the possibility that this variation may not be 
real, but an instrumental effect due to pileup. This is 
possible in principle because pileup increases with flux 
and serves to lower soft count rates and increase hard 
count rates, any residual pileup for which we have not 
accounted could cause a "hardening" of the spectrum as 
flux increases. However after detailed testing we con- 
clude that the variation is not due to pileup effects, be- 
cause: 

1. Modeling of pileup i n the Chand ra data using the 
ISIS pileup kernel ijPavigl 12001]) shows that the 
~50% changes in the flux that we observe can only 
change the power law slope by 0.01, which is much 
smaller than the observed variation in the Chandra 
spectrum. 

2. Variation in P is also observed in the XMM observa- 
tions, for which piled-up events have been removed. 
If we use all the events (including piled-up ones) in 
the X201 analysis, we find that the average value 
of P is flatter by ~0.1 due to pileup effects, but the 
variation in P of ±0.5 remains the same. 

Therefore the pulse variation in the power law index ap- 
pears to be independent of any pileup. We treat this 
effect as real and consider physical interpretations. 

Power-law index variation with pulse phase has been 
observed in a number of X-ray pulsars. W e searched 
the sa mple of bright pulsars in Table 2 of iHickox et alJ 
(j2004|) for those with existing pulse-phase spectroscopy. 
We find that many sources show variations in P, and 
these tend to be correlated w ith the pulse inten sity. Some 
sources, such as LM C X-4 ( Woo et al. ' 19961 , Cen X- 3 
(|Burderi et al.ll2000D. Vela X-1 (»Kretschmar et al.ll997j) . 
and 4U 1626-67 l|Pravdo et alJ 11979^1 . become harder 
with increasing intensity, as we have observed for SMC 
X-1. However, Her X-1 '(Ramsav et al. 2002), V0332+53 
Hunger et al.l ll992^. and 4U 1538-52 ( ClarkeTaDIIIll 
become softer with inc reasing flux. XTE JOlll. 2-7317 
ijYokogawa et al.ll2000(l . RX J0059. 2-7138 ijKohno et alJ 
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Fig. 7. — Diagram of the twisted disk model, showing the relevant 
parameters. The disk is viewed at elevation O^hs = 0°- Inset: the 
same disk viewed at B^ba = 90° , showing the disk twist angle <j)tw ■ 

1200(1 , and GX 301-2 l|Endo et al.ll2000D show no varia- 
tion in r. 

Because of this variety of behaviors, it is difhcult to 
present a universal physical picture for the pulse-phase 
changes in the power law spectra of X-ray pulsars. In 
the standard picture of spectral formation in luminous 
pulsar beams, the accretion column is stopped by a ra- 
diatively driven shock. If photons in the brightest part 
of the beam are generated in strongest shock, they will 
have been scattered with larger Compton y parameter 
and so preferentially emerge at higher energies. This 
would make the center of the pulse harder and the out- 
side softer, as observed for SMC X-1. A mechanism for 
a softening with pulse flux is not immediately clear. We 
would like to stress, however, that variation in F with 
pulse phase is very common, and so must be addressed 
by successful models of X-ray pulsar spectral formation. 

5. REPROCESSING BY A WARPED DISK 

The soft component pulses obser ved are rough l y con- 
sistcnt with the picture argued by ' Hickox et al.l l)2004f) 
and Ncilscn et al. ( 2004) that double-peaked power law 
pulses come from near the neutron star's accreting polar 
caps, while the soft pulses come from a much larger ra- 
dius, at the irradiated inner region of the accretion disk. 

We seek to test this possibility by comparing the ob- 
served hard and soft pulse profiles to a relatively simple 
model, in which a warped, twisted inner accretion disk 
is illuminated by a rotating X-ray pulsar beam. 

5.1. Illuminated disk model 

Our model does not describe the entire disk, but only 
the region close to the magnetosphere where the bright 
reprocessed component must be emitted. The disk model 
consists of a series of concentric circles with varying 
tilt and twist relative to each other with distance from 
the center (see Fig. 7. A similar shape has been in- 
ferred for the large-scale structure of the disk in Her 
X-1, from changes in the pulse profiles with superor- 
bital phase ((Hcott et al. 2000; Lcahv 2002). Warped and 
twisted structures have als o resulted from numeri cal sim- 
ulations of warped disks iWiiers fc Pringlelll99S ) and of 
gas- magnetosphere interaction ijRom'anoTOer^ 12003(1 . 



In our calculations we place this surface (we will refer 
to it simply as the "disk" ) at roughly the location of the 
magnetosphere (~10^ cm), and test if precession of such 
a surface around the disk axis can create changing soft- 
component profiles similar to those observed. The shape 
of the surface is determined by the radii and tilt angles 
of the inner and outer circle (fin, ^out, ^in and 0out), and 
the offset angle or "twist" between them (0tw)- Letting 
9d be the angle of the disk above the rotation plane at 
any point (r, 0), we use the simple approximation: 

V — T' 

Od{r, (t>) = -{Oin + (Oout - Oin) -^) sin (0 - 0tw), 

^out ^in 

which is valid for 9 < 30°. 

It is also necessary to define a beam pattern from the 
neutron star. Complex beam shapes have been proposed 
for some X-ray pulsars; for example Her X-l's emission 
has been model ed a s a central pencil beam surrounded 
by a fan beam (jBlum fc Kra us 2000), or even a reverse- 
pointin g fan beam, em itted from above the neutron star's 
surface (fScott et 31X 20001. However, since we are model- 
ing coarse pulse profiles with only 10 bins of pulse phase, 
it is impossible to constrain such complex models, so we 
use a simple beam shape. 

The beam is described by two two-dimensional gaus- 
sians, superposed on an isotropic component which is 
required since the hard pulse fraction for SMC X-1 is 
only ^50%. We define the beam pattern in spherical 
coordinates, with the poles aligned with the rotational 
axis and 6* = at the equator. We set the neutron star's 
rotation to be parallel to the disk axis. Each gaussian is 
specified by its width {ah), its angle from the rotational 
plane (Oh), and its azimuthal angle {(ph)- We also define 
a fan opening angle (0fan), which we set to zero to model 
"pencil" beams, and to > 45° to model "fan" beams. 

We set the peak intensity of each gaussian equal to 3 
times the intensity of the isotropic component. Note that 
this normalization is somewhat arbitrary; when we fit the 
output model profiles to the data we allow the intensity 
of the pulses to vary. The beam pattern is given by: 

Fhie, 0) = 1 + 3e-("i-«f-i)'/2'^? + 3e-(°2-«f-2)'/2<T,^ 

where ai^2 is the angular distance between {0,4>), and 
(^bi,2j 0bi,2)- This beam pattern is normalized to a total 
hard emission luminosity of Lx = 3 x 10'^^ ergs s^^. 

This beam pattern is then swept around the disk. We 
take the disk to be opaque and calculate the luminosity 
absorbed by each patch of the disk surface. We assume 
that all this energy is re-radiated as a blackbody spec- 
trum, with temperature 

T ={dL/{ATiadA)f/^. 

Only those regions of the disk that are directly visible 
from the neutron star are illuminated. Finally, we cal- 
culate which of the illuminated regions are visible to the 
observer, and are not blocked by parts of the disk in 
front. 

We assume that the heated disk immediately re- 
radiates its thermal energy, so that there is no time lag 
between illumination by the X-ray beam and emission 
from the disk. This assumption requires that both the 
light-crossing time for the system and the cooling time 
for the heated disk are much shorter than the pulse pe- 
riod. In our model the illuminated region is at i? ~ 10® 
cm, so the light crossing time is ^10 ms. 
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TABLE 2 
Parameters for disk models 



Model ranges Best fit 



Parameter 


Pencil beam 


Fan beam 


Pencil beam 


Fan beam 


Lx (ergs s ) 


o X iU 


Q vy 1 n38 
o X lU 


Q . , 1 n3S 
o X iU 


Q 1 a38 
O X iU 


i?in (10* m) 


0.8 


0.8 


0.8 


0.8 


Rout (10« m) 


1 


1 


1 


1 


Outer angle fiout (deg) 


30 


30 


30 


30 


Inner angle 0i„ (deg) 


10 


10 


10 


10 


Twist angle 0tw (deg) 


-90, 90, 135 


90 


90 


90 


fibl (deg) 


35, 45 


0, 45 


35 


45 


eb2 (deg) 


-10, 0, 10, 20 


0, 45 


-10, 10=" 





4>hi (deg) 








0, 





</'b2 (deg) 


210, 225, 240 


180, 210, 225 


225, 2W 


180, 210=- 


Beam half-width (deg) 


30 


15 


30 


15 


Fan opening angle Sf^n (deg) 





45, 60, 90 





60 


Observer elevation 9obs (deg) 


20 


20 


20 


20 



^These beam parameters are required to fit the X201 pulse profiles. 



A rough estimat e for the coohng timescale is given by 
et alJ l)200Cl() . and is simply the total thermal en- 
ergy of the heated disk divided by its luminosity. The 
disk's thermal energy is the volume of the heated region 
times its thermal energy density. If we treat the disk as 
a partial spherical shell presenting solid angle to the 
neutron star, we have 

where d is the penetration depth of hard X-rays into the 
heated region. We assume that most X-rays penetrate no 
more than a few Compton depths, so that nd 3/(Tt ~ 
10^5 cm-2. For Rbb = 10* cm, rj/47r - 0.1, and fcTeB = 
0.18 keV, we have Ebb ~ 10'^^ erg. The cooling timescale 
is simply ^cooi = -EeB/isoft, which for Lsoft ^ 10^"^ ergs 
s"i gives icooi ~ 10"^ s. 

Since both these timescales are much shorter than the 
pulse period, we can treat the reprocessing as instanta- 
neous. For each neutron star rotation phase, we calculate 
the total luminosity visible to the observer for the hard 
emission from the neutron star, appropriate to the ob- 
server's view of the given beam pattern. We also sum 
the observable reprocessed emission from the disk. We 
thus obtain "hard" and "soft" model pulse profiles that 
we can compare to observations. 

5.2. Model output 

We have created a number of illuminated disk mod- 
els, varying the disk shape (^in, ^out, 0tw), the beam pat- 
tern (6'bi, 6'b2, 0b2, 6'fan), and the observer's elevation an- 
gle ^obs- Because the model is computationally inten- 
sive, we did not perform a systematic study of the whole 
range of possible configurations. We initially ran the 
model for a broad variety of disk shapes and beam pat- 
terns, then limited our systematic study to those param- 
eters that could roughly reproduce the observed prop- 
erties of SMC X-1. These parameter ranges are given 
in Table 2. We set the outer disk to be inclined at 30°, 
which is consistent w ith dis k inclination estimates of 25°- 
58° iLutovinov ct all 120041) . and set the inner disk at a 
somewhat smaller inclination of 10°. This disk shape 
results in 20-30% of the neutron star luminosity being 
reprocessed. We set the elevation angle of the observer. 




Fig. 8. — Model twisted, tilted accretion disk as used in our 
reprocessing calculations. Shown is the model which allows the 
best fit to the observed (see Table 2). The shaded regions show 
areas on the disk which are illuminated by the central source and 
are visible to the observer. Lengths are in 10* cm. 



with respect to the neutron star's rotation plane, to be 
9ohs = 20°. If the neutron star spins in the plane of its 
orbit, this corresponds to an orbital inclination i = 70°. 
This is consistent with the estimated range of i ~ 65- 
70° from optical studies, assuming approximately Roche 
geometry for the coni panion star ijRevnolds et al.llT993t 
iHutchings et aljfl977|) . To obtain the observed black- 
body temperature of ~0.18 keV, we placed the reprocess- 
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ing region between rin — 0.8 x 10® cm and Tout = 1 x 10* 
cm. 

The model disks consisted of a grid on (r, </>) of 100 
elements in each dimension. For each disk shape we cal- 
culated the emission at 30 separate pulse phases and 8 
disk precession phases. The observer's views of the disk 
at several precession phases are shown in Fig. 8. In 
the upper left panel of Fig. 8, the neutron star has just 
emerged from behind the processing disk, which corre- 
sponds to the start of the high superorbital intensity 
state for which (j)so = 0. Therefore we set i/iproc = 
for this model disk orientation and define the other pre- 
cession phases accordingly. For each model beam and 
0prec: wc calculatcd the observed pulse profiles for the 
direct emission from the neutron star (the hard profile) 
and the reprocessed emission from the disk surface (the 
soft profile). 

5.3. Pencil beam models 

We begin by examining pencil beam shapes, for which 
^fan = 0. To reproduce the hard pulse profiles observed, 
we found that the two beams could not be anitipodal, 
that is 6»bi 7^ 0h2 and (/)b2 180°(with (^bi = 0°). The 
observed hard pulse peaks are not exactly 180° out of 
phase, so wc varied (/)b2 between 210° and 240°. To pro- 
duce the relative heights of the two hard pulses, we set 
6'bi to 35° or 45° and varied 6'b2 between -10° and 20°. 

While these represent only a small subset of the pos- 
sible disk/beam configurations, they allow us to test 
whether, in principle, reprocessing by a precessing disk 
can reproduce the observed pulse profiles. The model 
results for pencil beams show that: 

1. Most disk configurations display some kind of vari- 
ation in the shape and phase of the pulses in the 

soft, reprocessed component. 

2. The relative shapes of the hard and soft pulses vary 
most strongly with the shape of the illuminating 
beams, and less strongly with the detailed shape of 
the disk. 

For our array of models, we performed fits of the 
model pulses to the observed power law and blackbody 
flux profiles (shown in Fig. 5) for each observation. 
We allowed the normalization of the calculated pulses 
to vary, in order to account for any differences in pulse 
fraction between the observed pulses and our model. 

For each model we first fit the hard profiles to the 
observed power-law profiles, allowing the scale and the 
phase offset of the model profiles to vary. This achieved 
a reasonable fit (x^ ^ 2) for many pencil-beam models. 
We next attempted to fit the blackbody pulses by varying 
only the amplitude of the model profiles; the phase was 
fixed to that from the hard profile fit. For most model 
parameters, the shape and phase offset of the model soft 
pulses do not match those in the data, and no good fit is 
achieved. 

However, for one particular disk shape (Fig. 8, param- 
eters given in Table 2), we have achieved reasonable fits 
to the hard and soft profiles for all five observations (see 
Fig. 9). Between the observations we needed to change 
only: (1) the precession phase of the disk, and (2) for 
X201, the beam orientation (^b2 and (I)h2)- The reduced 



values for the individual fits are given in Table 3. 
Note that for the Chandra soft pulses, even the best fits 
have xt > ^ due to the complex shapes and small errors 
on the observed pulses. Qualitatively, however, the soft 
profiles have roughly the right features and phase offset, 
so we consider the fits satisfactory. The values in Table 
3 show that for each observation, there is a clear mini- 
mum in the values as a function of ipprec- This occurs 
because the reprocessed pulse profiles vary as the disk 
precesses and thus fit the observations for only specific 
disk orientations. 

The soft component pulses are fit best for c/iprcc = 
0.125 0.25 for Chandra, 0.375-0.5 for XlOl, and 0.75 for 
X201 (corresponding to the bold boxes in Fig. 9). For 
Chandra and XlOl these are consistent with the ^so es- 
timates of 0.16 and 0.37. There is not good agreement 
for X201, for which (pso = 0.42. However we note that 
X201 is also fit reasonably well for other phases includ- 
ing ^piec = 0.375, which would be consistent with its 
observed superorbital phase. 

5.4. Fan beam models 

We have also used our model to examine reprocessing 
of a fan-shaped pulsar beam, for ^fan > 45°. For the pen- 
cil beam we found that disk parameters do not strongly 
affect the reprocessed pulse profiles, so for the fan beam 
we fixed the disk parameters at the best-fit values for the 
pencil-beam case. We then varied the beam parameters 
with the values given in Table 2, setting the fan opening 
angles equal for both beams (^fani = ^fan2 = ^^fan)- In 
general, the fan beam is less successful than the pencil 
beam in reproducing the observed hard pulses, because 
in most cases the fan beam produces three or four peaks 
in the hard profile as both edges of each fan sweep past 
the observer. 

However, in the cases where we view the edge of each 
fan, or when there is only one fan ((?fan = 90°), we see 
only two hard pulses. For one beam configuration, we 
find fits that are roughly comparable to those for our 
best-fit pencil beam case (see Fig. 10 and Table 3). 
As for the pencil beam case, the fits to X201 require a 
slightly different beam configuration ((^b2 = 210° rather 
than 180°). 

Again we find that the best model fits occur for differ- 
ent ^prec for different observations, but in this case the 
best fits are for models with late disk precession phases 
((/-prcc = 0.625 for XlOl, 0.75 for Chandra, and 0.875 
for X201). This is not consistent with these observed 
4>so, which range from 0.16 to 0.42. In fact late preces- 
sion phases correspond to the low state of the superor- 
bital cycle, when the disk would be covering the neutron 
star (Fig. 8). Therefore we do not expect that these 
model parameters give a realistic description of the sys- 
tem. However, this result does suggest that in principle a 
fan beam, like a pencil beam, might produce the chang- 
ing soft profiles observed. 

6. DISCUSSION 

We have constructed a model of reprocessed pulses 
from a warped, precessing X-ray pulsar disk, and have 
succeeded in roughly reproducing the changing power law 
and blackbody pulse profiles observed in SMC X-1. This 
suggests that reprocessing by the accretion disk is a plau- 
sible explanation for the observed X-ray pulses in SMC 
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Fig. 9. — Hard and soft pulse profiles output by the disk reprocessing model, fitted to the observed profiles for the five observations. The 
disk and beam parameters are given in Table 2. Each column shows observed spectral component pulses (as in Fig. 5) from one of the 
five observations. The top row shows the power law profiles fitted by the pulses from the neutron star, and the lower rows show blackbody 
profiles fitted by the model reprocessed pulses. Reprocessing for 8 disk precession phases </<proc are shown. Note that the observed soft 
pulses vary in phase and shape with </)prcc- Highlighted in bold are the profiles that best fit for the observed power-law pulses (top) 
and reprocessed pulses (below). This disk and laeam configuration fits the Chandra data best at (jiprcc = 0.125—0.25 and XlOl best at 
<pprec — 0.375—0.5. The X201 pulses require a slightly different beam pattern (see Table 2) and fit the data well for (/>prcc = 0.625. Note 
that the X201 soft pulses are weaker than the other four, and are shown using a different scale. 



X-1, and is supported by the fact that other possible ori- 
gins of the soft component can be ruled out on physical 
grounds (Hickox et al. 2004). Because our model has a 
number of adjustable parameters, we cannot place real- 
istic constraints on the disk and beam geometries. How- 
ever, we can conclude that the long-term evolution of the 
SMC X-1 spectral behavior is consistent with precession 
of an illuminated inner disk. 

It remains unclear whether the long-term behavior of 
the power law and blackbody pulse profiles for SMC X-1 
is in fact correlated with the superor b ital p eriod. In an 
XMM study of Her X-1, iZane et all l)2004|) found that 
changes in the phase offset between the hard and soft 
pulses may vary monotonically with superorbital period, 
and so may be directly related to the disk precession. 
We are unable to say this for SMC X-1, because our de- 
termination of the superorbital phases is uncertain and 
because we only have data from a few epochs. Additional 
data would be very useful for solving this problem. Of the 
existing data taken since the launch of RXTE (which is 
required to determine ipso) there is only one high-state 
observation, with BeppoSAX, 2 March 1997 for which 
pulse phase-resolved spectroscopy has been performed 
i|Naik fc Padll2004a|) . This observation was taken at 
4>so ~ 0.5, and shows the familiar two-peaked power law 
profile and a wide blackbody pulse, offset from the power 



law pulses, which is broadly consistent with the observa- 
tions presented here. Other studies, such as pulse-phase 
spectroscopy performed with ASCA (Paulct al. 2002), 
show similar behavior but have no clear determination 
for ^so- 

A future observational program, comprising 6-8 expo- 
sures taken at regular intervals in 0so across the high 
state in a single superorbital cycle, would be ideal in 
allowing us to characterize the long-term evolution by 
studying gradual changes in the X-ray emission. In ad- 
dition, a few exposures at redundant (j)so in neighbor- 
ing superorbital cycles would help test the periodicity of 
these variations. It would be also interesting to perform 
a more limited study, similar to the present one, for LMC 
X-4. If that source shows long-term variation in the soft 
profiles which is also consistent with a disk reprocessing 
model, this will give support to a consistent picture of 
the accretion processes in this group of systems. 

For the parameter ranges we have chosen, the pen- 
cil beam model succeeds better than the fan beam in 
matching the pulse profile data at the known superor- 
bital phases. However the beam models are simple and 
there is much of parameter space that we have not ex- 
plored, which makes any preference of a pencil beam 
tentative at best. On theoretical grounds, the beam 
shape depends on the details of the accretion cohuim. 
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Fig. 10. — Same as Fig. 8, for a fan rather than a pencil beam shape. The profiles correspond to the fan beam configuration that best 
fit the observations, with parameters given in Table 2. 
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TABLE 3 

VALUES FOR MODEL PROFILE FITS^ 





C102 


C103 


C104 


XlOl 


X201 






Pencil Beam 








n Q 


0.8 


0.7 


o .u 


1 Q 


n nnn 
u.uuu 


1 ±.0 


15.9 


15.1 




1 A 


u. ±zo 


o .y 


10.2 


8.0 


O.U 








10.3 


7.9 


2 2 


-L .O 


U.O / 


Q 9 


13.4 


11.6 


1 




u.ouu 


10 7 
-LU. 1 


15.1 


14.0 




9 n 




1 n Q 


15.4 


14.1 


2 1 




0.750 


9.5 


13.9 


11.7 


2.2 


1.5 


0.875 


11.3 


15.9 


14.9 


2.6 


0.9 






Fan 


Beam 






Hardt- 


1.7 


1.2 


1.1 


0.7 


5.0 


0.000 


12.0 


16.9 


16.7 


4.0 


1.1 


0.125 


10.6 


15.4 


13.8 


3.5 


1.8 


0.250 


10.1 


15.2 


13.8 


3.1 


1.0 


0.375 


11.8 


16.5 


15.4 


3.6 


1.8 


0.500 


12.6 


17.6 


16.9 


3.9 


2.4 


0.625 


6.1 


7.3 


6.3 


0.6 


2.4 


0.750 


4.8 


6.2 


5.3 


1.4 


2.6 


0.875 


12.8 


17.5 


17.5 


4.0 


0.7 



^See Figs. 9 and 10 for model profiles. The hard and soft profile fits 
have 8 and 9 degrees of freedom, respectively. 



""For fits to the hard (power law) profiles. 



with polar cap emission from the neutron star's surface 
producing a pencil beam, and emission from the sides of 
the accretion co lumn pro ducin g a fan beam (e.g. Nagas3 
IT9X9tlMe.szar~ fc Riffe rtHTflS^ . Lx for SMC X-1 is close 
to or at the Eddington hmit (LEdd = 1.8 x 10^® ergs 
s~^ for a 1.4 Mq neutron star), which suggests that 
in the accretion column, where the accretion flow will 
be most intense, radiation cannot escape out th e top of 
the column and so a fan b eam is required (e.g. iBecked 
ll998HBecker fc Wolfal20n5j) . However, any reahstic pre- 
scription for pulsar beams involves a detailed descrip- 
tion of the accretion flow as well as effects such as grav- 
itational light-bending (Meszaros & Riffert 1988; Leahy 
1200.'^ ■ These can create complex beam patterns that 
may involve both pencil and fan shapes, and are beyond 
the scope of our simple model. 

While the present work cannot place strong constraints 
on the beam geometry, we note that our best-flt pulsar 
beam for both pencil and fan shapes is highly inclined 



with respect to the neutron star's rotation, and one beam 
is offset from the antipodal position of the other beam. 
This may suggest that the magnetic field is inclined 
from the rotation axis and is not exactly dipole. Dis- 
torted dipole fields ha ve also been sugges ted for H er X-1 
lIRIum Kra,u Fi [2nn?l and Cen X-3 (Kraus et a,].l il9961 
from detailed analysis of pulse shapes, and this possibil- 
ity will be interesting to test with future observations. 

7. SUMMARY 

We have performed pulse-phase spectroscopy of SMC 
X-1 from several epochs and have found that: 

1. The spectra are well fit by a soft blackbody plus 
hard power-law spectrum, and the soft and hard 
components vary separately with pulse phase. 

2. The shapes of pulses in the hard component remain 
essentially the same at different epochs, but the 
soft component pulses change on timescales > a 
few days. 

3. These varying soft profiles can be explained, in 
principle, by a model in which the soft emission 
comes from reprocessing of the rotating pulsar 
beam by a warped inner accretion disk that pre- 
cesses with time. Precession of such a disk is be- 
lieved to cause the known superorbital periodicity 
in SMC X-1. 

4. Further observations are required to establish 
whether the long-term changes in the hard and soft 
pulses of SMC X-1 are in fact correlated to the su- 
perorbital period, and to determine the precise disk 
and beam geometries. 
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